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ABSTRACT 

The effect of a non-trivial topology on the temperature correlations on the cosmic 
microwave background (CMB) in a small compact hyperbolic universe with volume 
comparable to the cube of the curvature radius is investigated. Because the bulk of 
large-angular CMB fluctuations is produced at the late epoch in low fio models, the 
effect of a long wavelength cut-off due to the periodic structure does not lead to the 
significant suppression of large-angular power as in compact flat models. The angular 
power spectra are consistent with the COBE data for fio^O.l. 
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1 INTRODUCTION 

Einstein's equations do not specify the global structure of 
spacetime. In other words, to a given local metric, a large 
number of topologically distinct models remain unspecified. 
In the absence of the unified theory that describes the global 
structure as well as the local one, one must resort to the ob- 
servational methods to determine the global topology of the 
universe. 

Assuming that the spatial hypersurface is homogeneous, 
the observed high degree of isotropy in the cosmic microwave 
background (CMB) points to the Friedmann-Robertson- 
Walker (FRW) models as the best candidates of the cos- 
mological models. However, if one would allow the spatial 
hypersurface being multiply-connected, a variety of locally 
FRW models which are globally anisotropic and inhomoge- 
neous may be consistent with the current observational data. 

Constraints on the topological identification scales us- 
ing the COBE data have been obtained for some flat models 
with no cosmological constant (Stevens, Scott & Silk 1993; 
de Oliveira, Smoot & Starobinsky 1996; Levin, Scannapieco 
& Silk 1998) and some limited compact hyperbolic (CH) 
models (Levin, Barrow, Bunn & Silk 1997; Bond, Pogosyan 
& Souradeep 1998). The large-angular temperature fluctua- 
tions discovered by the COBE constrain the possible num- 
ber of the copies of the fundamental domain inside the last 
scattering surface to less than ~8 for compact flat multiply- 
connected models. 

On the other hand, a large amount of CMB anisotropies 
on large scales could be produced in the low density universe 
due to the decay of gravitational potential near the present 
epoch (Cornish, Spergel & Starkman 1998). Therefore we 
expect that the constraint on the possible number of copies 



is less stringent for CH models. However, since the effect of 
the non-trivial topology becomes more and more significant 
as the volume of the space decreases, it is very important 
to investigate the viability of the CH models with small co- 
moving volume. 

From a theoretical point of view, the "smallness" of 
the spatial hypersurface is an advantage for giving a nat- 
ural mechanism leading to homogeneity and isotropy. It is 
well known that geodesic flows on CH spaces are strongly 
chaotic. Therefore, initial perturbations would be smoothed 
out due to the mixing effects (Lockhart, Misra & Prigogine 
1982; Gurzadyan & Kocharyan 1992; Ellis & Tavakol 1994). 
In inflationary scenarios, a certain physical process is in- 
dispensable that homogenises the initial patch beyond the 
horizon scale before the onset of inflation for accomplishing 
the sufficient smoothing of the observable universe (Gold- 
wirth & Piran 1989; Goldwirth 1991). The chaotic mixing 
in CH spaces may provide a solution to the pre-inflationary 
initial value problem (Cornish, Spergel & Starkman 1996). 

If we live in a small universe which is defined to be a 
locally homogeneous and isotropic space that is multiply- 
connected on scales comparable to or smaller than the hori- 
zon, the future astronomical satellite missions such as MAP 
and PLANCK might reveal some specific features in CMB 
(Cornish, Spergel, & Starkman 1998; Weeks 1998). 

So far, a variety of CH manifolds have been constructed 
by mathematicians. However, the number of the known CH 
manifolds with small volume is relatively small. In this pa- 
per, we investigate CH models whose spatial hypersurface 
is isometric to the Thurston manifold which is the second 
smallest in the known CH manifolds with volume 0.98139 
times cube of the curvature radius. The smallest one is 
the Weeks manifold with volume 96 percent of that of the 
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Thurston manifold(see e.g. Fomenko & Kunii 1997). How- 
ever, the fundamental domain(which tesselates the infinite 
space) of the Thurston manifold is much simpler than that 
of the Weeks manifold. For simplicity, we investigate the 
Thurston models rather than the Weeks models. The fun- 
damental domain of the Thurston manifold is a polygon with 
16 faces, which can be constructed by appropriately identi- 
fying 8 faces with the remaining 8 faces(see the appendix 
of Inoue 1999a). It should be noted that the volume of CH 
manifolds must be larger than 0.16668 times cube of the 
curvature radius although no concrete examples of mani- 
folds with such small volumes are known (Gabai, Meyerhoff 
& Thurston 1996). 
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Figure 1. Number function and the Weyl's asymptotic formula. 



2 COMPUTATION OF EIGENMODES 

So far various kinds of numerical techniques have been pro- 
posed to overcome the difficulty of computing the CMB 
in CH models. For several CH models, CMB fluctuations 
have been computed using the method of images with- 
out carrying out the mode expansion (Bond, Pogosyan, & 
Souradeep 1998). They obtained the result that the COBE 
data strongly constrains the CH models so that the comov- 
ing volume of the fundamental domain are at least com- 
parable to the comoving volume inside the last scattering 
surface. Since the method of images requires the sum of ex- 
ponentially increasing images, it is difficult to obtain the dis- 
tinct eigenmodes which are necessary to estimate the effect 
of the power spectrum with discrete peaks. Alternatively, 
one of the author proposed a numerical approach called the 
direct boundary element method for computing eigenmodes 
of the Laplace-Beltrami operator (Inoue 1999a). 14 eigen- 
modes have been computed for the Thurston manifold. It is 
numerically found that the expansion coefficients behave as 
if they are random Gaussian numbers. 

In this work, we have numerically computed 36 eigen- 
modes in the Thurston manifold up to k — 13 (the curva- 
ture radius is normalized to one) which are approximated by 
quadrature shape functions which converges to the solutions 
faster than constant valued shape functions. As we shall see, 
the contribution of the higher modes to the angular power 
spectra on large angular scales are relatively small for low- 
density models. In other words, the effect of the non-trivial 
topolgy is almost determined by the lower modes. We con- 
firm the previous computed eigenvalues within \Sk\ ^ 0.01. 

We see from figure 1 that the number of eigenmodes 
below k is nicely fitted to the Weyl's asymptotic formula 
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Vol(M)(fc 2 
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where Vol(M) denotes the volume of a manifold M. The 
random Gaussian behavior is again observed for 31 modes 
5.404 ^ k < 13 but five degenerated states have an eigenmode 
which shows the non- Gaussian behavior due to the global 
symmetry of the fundamental domain. It is found that the 
five eigenmodes have Z2 symmetry (invariant with respect 
to the rotation by an angle n) on the center (where the min- 
imum length of the periodic geodesic which lies on the point 
is locally maximal) of the fundamental domain. In this case, 
one would observe an axis around which the fluctuation is 
rotationally symmetric at the center. Therefore, the correla- 



tion between expansion coefficients leads to a non-Gaussian 
behavior. Nevertheless, it is found that appropriate choices 
of the linear combination of the degenerated modes recover 
the generic Gaussian behavior. Furthermore, the symmetry 
of CH manifolds depends on the observing point. If one ran- 
domly choose a point on the manifold, the probability of 
observing an exact symmetry of the manifold is very small. 
The result supports the previous investigations of the ex- 
pansion coefficients which show the Gaussian behavior in 
classically chaotic systems (Aurich & Steiner 1989; Haake & 
Zyczkowski 1990) although the global symmetry in the sys- 
tem can hide the generic property (Balazs & Voros 1986). 



3 TEMPERATURE FLUCTUATIONS 

Perturbations in CH models can be written in terms of linear 
combination of eigenmodes on the universal covering space 
multiplied by the expansion coefficients and the initial fluc- 
tuations plus time evolution of the perturbations. The ex- 
pansion coefficients include the information of the periodic- 
ity in the universal covering space. As CH models are locally 
homogeneous and isotropic, the time evolution of the per- 
turbations coincides with that in open models. 

The dominant physical effects producing CMB 
anisotropies (Hu, Sugiyama & Silk 1997) on large angular 
scales are the ordinary Sachs- Wolfe (OSW) effect (Sachs & 
Wolfe 1967), which is the gravitational redshift effect in be- 
tween the last scattering surface and the present epoch, and 
the integrated Sachs- Wolfe (ISW) effect, which is the grav- 
itational blue-shift effect caused by the decay of gravita- 
tional potential at the curvature domination epoch, l + 2~ 
(1 — flo)/flo- For the COBE scales, we can ignore the contri- 
bution from the acoustic oscillations. Then the time evolu- 
tion of the adiabatic growing mode of the Newtonian grav- 
itational potential is analytically given as (see e.g. Kodama 
& Sasaki 1986; Mukhanov, Feldman & Brandenberger 1992) 



$ t (77)=$ t (0) 



5(sinh 2 77 — 3?7 sinh rj + 4 cosh T7 — 4) 
(cosh 77 — l) 3 



(2) 



where r\ denotes the conformal time. The two-point tem- 
perature correlations in a CH cosmological model can be 
written in terms of the gravitational potential. Assuming 
that the initial fluctuations obey the Gaussian statistic, and 
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neglecting the tensor-type perturbations, the angular power 
spectrum Ci can be written as 



(2l + l)Ci = ^2(\a lm \ 2 ) 



E; 



4tt 4 V*{v) 
{v 2 + l)Vol(M) 
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where 
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Here, v = \/k? — 1, "P^iv) is the initial power spectrum, and 
77* and n are the conformal time of the last scattering and 
the present conformal time, respectively. X„i denotes the 
radial eigenfunctions in open models and £wm denotes the 
expansion coefficients. From now on we assume that the ini- 
tial power spectrum is the (extended) Harrison-Zeldovich 
spectrum i.e., V<s>(y) = Const.. 

Although the low-lying modes give an appreciable con- 
tribution to the large angular power, contributions of higher 
eigenmodes may not completely be negligible. While the 
computation of highly-excited eigenmodes is a difficult task, 
we have so far succeeded to calculate the exact eigenmodes 
up to k — 13 as we mentioned before. However, we are go- 
ing to assume that £„; m 's are also random Gaussian num- 
bers for higher modes. Since the information of the peri- 
odicity in the real space is lost by this approximation, we 
will only employ this approximation to the statistics in the 
fc-space which is expected to be not changed because the 
periodicity is not apparent in the fc-space. As CH models 
are globally inhomogeneous, the expected correlation statis- 
tics depend on the point of the observer. Therefore, one can 
interpret that one realization for the expansion coefficients 
corresponds to a certain point of the observer in the fun- 
damental domain. In order to apply the random Gaussian 
approximation, one must also estimate the variance of the 
expansion coefficients. The expansion coefficients are writ- 
ten in terms of eigenmodes u v and spherical harmonics Y; m 



(5) 



ivi m X vl (xo) - \ u„(xo,0,</>)Y lm (6,<l>)dtt. 



It should be noted that (pi) is satisfied at arbitrary ra- 
dius Xo- Let us consider a sphere with large radius \o > 
> 1 on the Poincare ball which is the image of the up- 
per hyperboloid in the four-dimensional Minkowski space 
(1/0,2/1,3/2,2/3) by a stereographic projection onto the unit 
ball on the (0,2/1,2/2,2/3) plane using a point (—1,0,0,0) as 
the base point. One can expect the random behavior of the 
mode functions on the sphere as the surface of the sphere 
which is pulled back by the discrete isometry group fills 
the fundamental domain ergodically. The (apparent) angu- 
lar fluctuation scale 59 of fc-mode is approximated in terms 
of two parameters \o and k as , 



60- 



16tt 2 Vol(M) 



fc 2 (sinh(2(xo + r ave )) - sinh(2(xo - r ave )) - Ar ave ) 



(6) 



where r ave denotes the averaged radius of the inradius and 
outradius of the fundamental domain. One can approximate 
u v'iXo) ~ u v{x'o) by choosing an appropriate radius Xo 
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Figure 2. ST t /T = y/l{l + 1)C 1 /{2tv) for the Thurston models 
Qq = 0.2 (left) and Co = 0.4 (right) using expansion coefficients 
derived from 36 eigenmodes only(stars) and that using these co- 
efficients and random Gaussian numbers for 13 < k < 50 with 100 
realizations (diamonds). Eigenvalues for higher modes arc approx- 
imated by Weyl's asymptotic formula. 




Figure 3. A simulated sky map of the microwave background 
(convolved with the COBE DMR beam) in the Thurston model 

n = 0.2. 



which satisfies k 2 exp(— 2x' ) ~ k' 2 exp(— 2\o). Averaging 
(n) over I and m, one obtains 



{%■ 



,2\ exp(-2xo) /| C 
l| /~cxp(-2xo)\ IW 



2 >- 



(7) 



". We have found that the com- 



which gives (|£„im| 2 ) ~ 

puted variances of £„; m 's for 2 sj I ^ 20, —l^m^l are 

remarkably in good agreement with the analytical estimate. 

From figure 2, one can see that the uncertainty in the 
Gaussian approximation is very small. Remarkably, each re- 
alization gives almost the same value so that 100 points for 
given I are plotted as a tiny speck. The contribution of higher 
modes becomes significant as flo is increased because the 
curvature dominant era is shifted to the late time so that 
the OSW effect becomes dominant over the ISW effect. It 
is found that contributions of the modes k > 13 to C; for 

2 ^C I ^ 20 are approximately 7 percent and 10 percent for 
fio — 0.2 and f2o = 0.4, respectively. Thus contribution of 
modes k > 13 which we employ Gaussian approximation is 
almost negligible on large angular scales especially in low fio 
models. One realization (for the initial fluctuation) of a typ- 
ical CMB fluctuation as seen by COBE is plotted in figure 

3 for fio = 0.2. In the simulation, we used only "exact" 36 
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Figure 4. 5T t /T = <Jl{l+ l)Ci/(2ir) for the Thurston mod- 
els at center (diamonds) and ensemble averaged values (boxes) 
and open models (stars) with Ho = 0.2 (left) and Qq = 0.4 
(right). Because of the global inhomogeneity, STi/T's have de- 
pendence on the observing points inside the fundamental domain 
that causes the uncertainty in STi/T. Two-sigma "geometric vari- 
ance" is shown in vertical lines which has been obtained by 500 
realizations for the expansion coefficients. 
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eigenmodes. We have chosen a point where the injective ra- 
dius is maximal as the center (belonging to the "thick" part 
of the manifold). One can see that the structure due to the 
periodical boundary conditions is not apparent. However, 
approximated number of copies of the fundamental domain 
inside the last scattering surface is ~ 500 for the Thurston 
model with fio =0.2. Therefore, the effect of the non-trivial 
topology is expected to be significant. 

The mode cut-off at k — 5.404 which corresponds to the 
largest wavelength inside the fundamental domain causes 
the suppression of the angular power on large angular scales 
as in compact flat models. However, the decay of the New- 
tonian potential in the curvature dominant era makes the 
difference. Since the bulk of the large angular power comes 
from the decay of the potential well after the last scattering 
time, the large angular power does not suffer the significant 
suppression. We see from figure 4 that the slope of the large 
angular power is not steep even for the model with Qq = 0.2 
in contrast to the compact flat models without cosmological 
constant. The two peaks in the power spectrum for the CH 
model are important in understanding the effect of the non- 
trivial topology. The angular scale which gives the first peak 
is equivalent to the angular fluctuation scale of the lowest 
eigenmode (A; = 5.404) on the last scattering surface. Substi- 
tuting the comoving radius of the last scattering surface in 
unit of the curvature radius R cur v, 



Rlss = Rcurv cosh 1 (2/Q - 1) 



(8) 



into (n) gives the angular scales I = 17 for flo = 0.2 and 
I = 7.4 for fi = 0.4. Beyond this scale, the OSW contribu- 
tion is strongly suppressed as in compact flat models. How- 
ever, eigenmodes with angular scales below the given scale 
at the last scattering can have large angular scales after the 
last scattering. Therefore, in the presence of the ISW ef- 
fect, the suppression of the power beyond the scale which 
corresponds to the first peak is very weak in contrast to 
flat models. The angular scale which gives the second peak 
corresponds to the scale of the projected lowest eigenmode 
at the last scattering. Below this scale, the angular power 



Figure 5. STi/T = y/l(l + l)Ci/(2n) in fjK for the Thurston 
model with f2o = 0.1 (top left) and Cq = 0.4 (top right) and 
Q = 0.2 (bottom left) and fl = 0.6 (bottom right) The light- 
gray band corresponds to the one-sigma cosmic variance in center 
and the darkgray band corresponds to the net variance of the 
one-sigma cosmic and the one-sigma "geometric" variances. The 
COBE DMR measurements analysed by Gorski and Tegmark are 
plotted in diamonds and boxes respectively. 



asymptotically converges to that of open models because the 
effect of the modes with wavelength larger than the cut-off 
wavelength is negligible. Since we have ignored the effects 
of subhorizon perturbations at the last scattering such as 
the so-called 'early' ISW effect during the matter-radiation 
equality epoch and the Doppler effect due to the acoustic 
velocity, the angular power on large to intermediate scales 
must be slightly boosted. However, these effects are irrele- 
vant to the global effect of the non-trivial topology inasmuch 
as one considers the typical topological identification scale 
that is not significantly smaller than the present horizon. 

In figure 5, the angular power spectra for low-Omega 
models are plotted with the COBE data (Gorski et al, 1996) 
(diamonds) . They have been calculated using 36 eigenmodes 
and the Gaussian approximation taking account of ~10 per- 
cent contributions from higher eigenmodes. The slope of the 
power becomes steep as fio is lowered since the ISW contri- 
bution transfers to the large scales. 

We have performed a simple \ 2 fitting analysis to the 
COBE DMR band power measurements (Tegmark 1997) 
(boxes) which are uncorrelated . We have adjusted the nor- 
malization of the initial power to minimise the value of x 2 ■ 
As shown in table 1, the angular power for a model with 
f2 = 0.1 is still within the acceptable range. The apparent pri- 
mordial spectral index is approximately n— 1.6 for f7 = 0.1. 
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n 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


x 2 
Q 


10.6 
0.15 


6.29 
0.51 


6.42 
0.49 


4.21 
0.76 


4.33 
0.74 


5.20 
0.64 



Table 1. \ 2 an d the probability Q that x 2 should exceed a par- 
ticular value by chance assuming that x 2 obeys the chi-square 
distribution with 7 degrees of freedom. 



4 CONCLUSIONS 

Thus the Thurston models with fi^O.l are not constrained 
by the angular power spectrum from the COBE data, which 
confirms the preliminary result by one of the author (In- 
oue 1999b). The peak at I ~ 4 in the COBE data may be 
merely the coincidence due to the large cosmic variance but 
it is interesting that a model with fio ~ 0.6 has the first 
peak in this scale. Consequently, the Thurston models agree 
well with the COBE data than any FRW models. The sim- 
ilar conclusion that the constraints Q ^ 0.3 for an orbifold 
model with volume 0.7173068i?c Url , have been obtained in 
(Aurich 1999). Although orbifolds have singular points, the 
behavior of eigenmodes for orbifolds is expected to be simi- 
lar to that of manifolds. Therefore, the result for an orbifold 
model supports our conclusion. 
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